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^ : 1 Introduction. 

The method of Supersymmetric (SUSY) Quantum Mechanics [Ij provides a very effective 
approach for investigation of different problems in Quantum Mechanics. The SUSY inter- 
twining relations [2] and the notion of shape invariance [3] are most powerful tools. In 

^ ' one-dimensional SUSY QM they jointly reproduced in a very elegant way all exactly solv- 
ed '. 

able models which were known by that time [1], [5]. In addition, the intertwining relations 
provide the wide class of the so-called quasi-exactly-solvable (QES) or, equivalently, par- 
tially solvable models |6]. Shape invariance for models invariant under a centrally extended 
superalgebra with an additional symmetry generator was demonstrated in 0. The shape 
invariance in the context of Calogero-like A^— body systems was studied in [8]. 
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In two-dimensional generalization of SUSY QM, the intertwining relations [9], [TO], [TT] 
and shape invariance [12], [13], [H], [15] also play very important role, but their realization 
is slightly different. First of all, it must be taken into account that different two-dimensional 
generalizations of SUSY QM exist. The direct generalization [9], [lO] with first order super- 
charges intertwines Hamiltonians with different matrix dimensionality. We shall consider 
another opportunity for the two-dimensional case - intertwining of two scalar Hamiltonians 
H^^\ ^(2) ]-,y supercharges of second order in derivatives [H], [13] : 

H('\x)Q+{x) = Q+{x)H^^\x), (1) 
Q-(f)if(i)(x) = H^^\x)Q~{x), (2) 

i7(l).(2) = _ Q2 ^ y(l),(2)(^). (xi,X2), 

where supercharges are chosen with Lorentz metrics in second order derivatives: 

Q+ = Ad+d. + 4C+9_ + AC^d+ + B{x)- Q- = {Q+Y; x± = Xi± X2. (3) 

As was shown in [TT], this choice of the metrics simplifies essentially the solution of inter- 
twining relations (111), ([2]). In the case of arbitrary form of second order part in Q"^, these 
relations are equivalent to the very complicated system of nonlinear partial differential equa- 
tions for the functions V"*^^^'^^-*, C±, B. For the case of Lorentz form ([3]) of second order part, 
this system of equations leads to C± = C±{x±) and to the expressions for the potentials (we 
remark that coefficient functions C± differ from that in previous papers [TT], [13], [12] by a 
multiplier 4): 

= 2(C; + Cl + C'_ + C2 ) + f,{x2) - = 
^ v^^\x+) + v^}\x_) + Mx2) - (4) 
= 2iCl -C'^ + Cl- C'_) + h{x2) - ^ 
vf{x^) + v^'\x^) + h{x2) - hix,), (5) 
B{x) = 4C+(x+)C_(x_) + + /2(X2), (6) 

where the functions C±{x±) and F{x) = fi{xi)+f2{x2) have to satisfy the following equation: 

9_(C_F) = -9+(C+F), (7) 
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the last (unsolved yet) equation of the above mentioned system of equations. This equation 
was solved in [11], [13] by means of choosing some simplifying ansatzes. Some of these 
solutions obey [H], [15] two-dimensional shape invariance property. 

In the present paper we shall study the general question: is there any, wide enough, 
additional class of shape invariant systems? Therefore, it is timely to remind for the reader's 
convenience the idea of shape invariance in SUSY Quantum Mechanics. 

In one- dimensional Quantum Mechanics with first order differential operators Q^, one 
refers to shape invariance [3], [1], when the partner Hamiltonians in ([1]), ([2]) both depend on 
some (multi) parameter a, and they have the similar shape, i.e. they satisfy: 



where a = a{a) is some new value of parameter, which depends on a, and TZ{a) is a (c- 
number) function of a. In the case of absence of spontaneous breaking of supersymmetry, 
this property allows to construct the analytical expressions for all wave functions \l/(x,a) in 
a pure algebraic way (see details in [Ij). The whole energy spectrum of the Hamiltonian is 
found as well. Among several relations of SUSY algebra intertwining relations ([T]), ([2]) play 
a crucial role in this approach. 

The notion of shape invariance was generalized onto two-dimensional SUSY Quantum 
Mechanics in [12], [13], [H], [15]. In this case, the partner Hamiltonians depend also on 
(multi)parameter a, the main relation ([8]) has the same form, but operators are of second 
order in derivatives. Thereby, in contrast to one- dimensional situation, have many zero 
modes, and one has chance to find only a part of the spectrum and the corresponding wave 
functions (quasi-exact-solvability) [12], [13]. This does not prevent us from finding the whole 
spectrum (exact solvability), using however some additional tools [T6] . 

The structure of the paper is the following. The general conditions which guarantee 
the shape invariance of two-dimensional potentials are in Section 2, using their integrability 
and the results on shape invariant one-dimensional systems. The possible ansatzes for the 
coefficient functions of the supercharges are defined. In Section 3 the explicit expressions 
for potentials are derived. The particular forms of these potentials for specific values of 
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parameters are given in Appendix. Some of potentials are recognized as already known, but 
others are new to the best of our knowledge. 

2 Two-dimensional shape invariance. 

The problem of constructing of most general form of shape invariant potential seems to be a 
rather difficult task. Even in one-dimensional SUSY Quantum Mechanics it is not yet fully 
solved (see discussion in [1], [5]). Only a class of such potentials was built, and it coincides 
with variety of well known exactly solvable potentials. It was proven recently [5] that no 
additional shape invariant potentials with so called additive shape invariance exist, outside 
this class. 

It is clear that one has not too many chances to solve the analogous problem in two- 
dimensional situation. Therefore, we do not pretend to find all existing shape invariant 
two-dimensional potentials. At best, we can find a wide class of such new potentials. It 
is necessary to remind that some shape invariant potentials were already obtained [Hj, 
[13] in the framework of polynomial SUSY QM. Among these potentials two-dimensional 
generalizations of Morse potential and Poschl- Teller potential must be mentioned specially 

m- 

It is convenient, without loss of generality, to choose parameters in such a way that the 
shape invariance condition ([8]) links Hamiltonians H^^''^\x]a) with difference (step) between 
parameters equal to unity, i.e. 

H^^\a + l;x) - H^^\a;x) = const. (9) 

It is evident from the explicit expressions (jlj), ([5]) for potentials, that the shape invariance ([9]) 
is equivalent to a pair of one-dimensional shape invariance with "superpotentials" C±{x±) : 

w£V± + 1; - v±\a±; x±) = c±, (10) 

i.e. Hamiltonians H^^''^^ obey actually two-parametric shape invariance with the two- 
component parameter a = (a+, a_). Besides the potentials V-^'^^(a;x) contain also 
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the terms f2{,X2) — which prevent separation of variables. 

The dependence of H^^''^^ on a seems superficially to be a generic dependence on two 
independent variables a+, a_. However, as will be clear from the following Eq. lfT^ there is 
a constraint. It is well known [TT], [13], that arbitrary Hamiltonians H^^\a; x), H^'^\a; x), 
which participate in SUSY intertwining relation, are integrable, i.e. they obey the symmetry 
operators (here of fourth order in derivatives): 

[R'''^\a;x),H^^\a;x)] =0, R'^'^'' {a; x) = Q+ {a; x)Q- {a; x) 
[R^'^\a;x),H^^\a;x)] = 0; R^^\a;x) = Q-{a;x)Q+{a;x). 

The shape invariance relation (Q allows to conclude that: 

[R^^\a + l;x)-R^^\a;x),H^^\a;x)] =0. (11) 

Although the operator R^^\a + l;x) (a; x) seems to be a symmetry operator for 

(a; x), one can check that it is of second order in derivatives. The symmetry operator of 
second order in momenta signals that the Hamiltonian is amenable to separation of variables 
[18] . [TT] . Since from the very beginning, we do not consider systems with separation of 
variables, we shall concern ourselves only with the case when this operator is a function of 
H^'^\a-x), i.e.: 

(a + 1; f ) - (a; x) = fxH^^'> (a; x) + u, (12) 

with fi, I/- const ants. 

Using expressions for Q"^, the l.h.s. of (fT2l) can be calculated explicitly: 

R^'\a + 1) - = ^s(^yW^a+ + 1) - vf{a+)^d^_ -8(^v^'\a_ + 1) - v^^\a^)^dl + 

+16( C;(a+ + l)C'_{a_ + 1) - C'{a+)C'_{a_)]+A( C+{a+ + l)d^B{a + 1) + 



+C_(a_ + l)d+B{a + 1) + C+{a+)d-B{a) + C_(a_)9+5(a) +B\a + 1) - B\a), (13) 



where the dependence on space coordinates was omitted for brevity. Comparing with (fT2|) . 
we obtain the necessary condition, which constants c± in the r.h.s of Eq. (fTO|) have to satisfy. 



Namely they should coincide: 

v^^\a+ + 1) - (a+) = t;i'^(a_ + 1) - v'i'\a^) = c. (14) 

According to analysis of [5] , where all known one- dimensional shape invariant potentials 
were reproduced, we shall consider the following form of dependence of " superpotentials" 
C± on parameters a± : 

C±ix±) = a±p±{x±) + r±{x±) + q±{a±), (15) 

where a±. are parameters of shape invariance with a unit step a± = aj. + 1, functions 
P±{x±),^±{x±) do not depend on a±, and q±{a±) do not depend on x±. It was shown in 
[5], that it is sufficient to use the two-term variants of ( fTSi) . therefore we shall consider 
below separately two different kinds of possible dependence of C±(x±) on shape invariance 
parameters: 

J) C±(a±) = a±p± + r±; (16) 

II) C±(a±) = a±p± + g±(a±). (17) 

We shall explore these ansatzes together with additional requirements, which are specific for 
two-dimensional systems, namely Eq.([7]) and f|T2l) . 

3 Construction of shape invariant potentials. 

The shape invariance of first kind can be considered, starting from Eq.f ll4p . which for the 
choice f|T6|) reads as: 

(2a± + 1)(p1 +p'±) + 2p±r± + 2r'^ = c/2, (18) 

where c is an arbitrary constant, the same for signs ± in the l.h.s. It follows that in this 
case, p±{x±), r±{x±) satisfy the system of differential equations: 

2,1 \ 2 I / J 

P±+P± = A±; p±r± + r^ = d±, 
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and its general solution is: 



P± = — 7 — r; 2;±(a:±) = a±exp (A±a;±) + 5±exp (-A±x±); (19) 
z±{x±) 

= ;^(«± + ^±/'-±(^±)^^±)=;^(«± + |4(^±)), (20) 

with integration constants a±. The form ( |T6l) of C± and expressions ( |T9l) . ( l20l) in terms of 
give an opportunity to take = Q hy means of transformation of parameter to 
(a-i- + (i±/A^). Thus, we shall continue with d± = Q below. Also, the direct calculations give 
that the constant c in the r.h.s. of Eq.( fT4l) is: c = 2(2a + 1)A^, and therefore A+ = A_ = A. 

Now we can go to the system of equations (JTj). Since the function F does not depend on 
a±, it follows from ([7j) for independent parameters a_|_ and a_, that the following relations 
must be fulfilled: 

d+{Fp+) = d-{Fp^) = 0; 9+(Fr+) + (9_(Fr_) = 0. (21) 

From the first two equalities, one can conclude that F = Const/ (p+(x+)j9_(x_)). This is just 
a case of factorizable function F, which was studied earlier in [Tl] in a general form. One 
of two solutions, which were found there, indeed obeys shape invariance, but it is amenable 
to standard separation of variables. By this reason, we shall not consider this case further, 
restricting ourselves to a case with one independent parameter instead of two: a_ = a_|_ = a. 
As a consequence, the possible shape invariance is still reduced to a pair of one-dimensional 
shape invariance in variables x±, but the restrictions onto functions become much less strong 
than ( 12T|) now: 

a+(Fp+) + 9_(Fp_) = 0; (22) 
9+(Fr+) + 9_(Fr_) = 0. (23) 

If one of the constants a+ or a„ in ( l20l) vanishes, it is clear from ( l23l) that F is factorizable, 
therefore we shall study two other choices of constants: 

la) = a+ = 0; lb) a^a^ ^ 0. 
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la). For this option, Eq. (!23|) is satisfied identically, but we have no direct way to solve 
Eq.( l22|) in a general form. We shall act in an indirect way. It follows from ( !T3|l . that 

R^^\a + 1) - = 8\'^{2a + 1) (^H{a) + 2\'^{2a + 1)^ + 

+2(2a + 1) (^4A2(/i - f,) + (p+ + + {p_ - p^)f, + 4p+p_(/i + f,)^ + 

+4 f (r+ + r_)/( + (r_ - r+)/^ + 2(p+r_ + p_r+)(/i + /s)^ . (24) 



According to Eg. (1121) and since the functions /i,2,P± do not depend on parameter a, each of 
last two terms in ( 12^ must be constant. In particular, 

4A2(/i - /s) + (p+ + p-)/( + {p- - p+)f2 + ^P+P-ifi + /2) = 2a;; a; = const. 

Together with (l22l) . this equation gives: 

4AVi(a;i) -u; + (p+ +p-)/{(xi) 

and after differentiation over xi, we obtain the equation for the function 



/;' + 5^^^i^/! + 8AV. = 2u,. (26) 

P++P- 



From definition (IT^ of p± : 

6(A^+p+p,) _ 3^ 

P+ +P- Zi'' 

where zi is defined as: 



(27) 



2;i(xi) = (T+(T_ exp (2Axi) — 5+5- exp (— 2Aa;i); 2;^' = AX^zi. (28) 
Eq.(l271) transforms (1261) as follows: 



whose solution is: 

oj kiz[ + ^2 
4A2 ^ 7 



/i(-i) = l^ + ^^^^, (29) 



with ki, k2— constants. Substitution back into (125|) gives: 

oj , kiz'2 — k2 



/2(x2) = -^ + ^^|-^, (30) 

with the definition of Z2 '■ 

Z2{.X2) = cr+(5_ exp (2AX2) — (5+cr_ exp (— 2AX2); Z2 = 4:X^Z2. (31) 

It follows from (!29l) and (130|) . that we may take the value a; = 0, since /i and /2 are defined 
up to an additive constant with opposite sign. 

Thus, the option la) gives the following shape invariant potentials: 

a-6- \ A;i((j+(T_ exp(2Aa:;i) + exp(— 2Axi)) + k2 

(cr_ exp(Aa;_) + 5_ exp(— Aa;_))^/ (o"+o"- exp(2Axi) — exp(— 2Axi))^ 
^ fci(o-+(5- exp(2Ax2) + exp(-2Ax2)) - A:2 
(a"+(5_ exp(2Ax2) — (5+cr_ exp(— 2Ax2))^ 

The particularized expressions for various choices of the arbitrary constants are given in the 
Appendix. 

lb). In this case we shall solve the system of equations ( 122|) . (!23|l directly. This task is 
simplified by assuming a+ = a_ = a in (120|) without loss of generality. This is possible due 
to homogeneity of equations f l22p . f l23p under multiplication of r± by constant factor. After 
simple manipulations we obtain the system (1221) . ( I23l) in the form: 



r'_p_ - p'_r_ + r+p'+ - r^p+ A2(r+ - r_) 

ai(lnF(p+r_ — r+p_)) = — = 

]9_i_r_ — r^p_ p^r_ — r^p^ 

r'_p_ - p'_r_ - r+p'_^ + r'^p+ X^{r-+ + r_) 

a2(lnF(j9_|_r_ — r^p_)) = — = 

p_i_r_ — r_|_p_ p_|_r_ — r+p_ 

and using expressions (fT9|) . ( 120!) : 

ai(lnF(p+r_ - r+p_)) = -A^^f^ = -^iHl!^ = -^^i ln(y; - y'_). 

y+ -y- y+- y- 

a2(lnF(p+r_ -r+p_)) = = -^±jtl_ = ln(y+ - y^). 

2/+ - 2/- 2/+ - 2/- 



These equations can be integrated explicitly: 

17 _ y+y- 



(33) 



{y'+-y'-r 

but it is necessary to take into account additionally that F = fi{xi) + f2{x2)- This gives 
restriction onto parameters in f|T9|) : cr+(5+ = cr_5_. Then, for 5+ 7^ 0, we obtain from fl33|) . 
that: 



/: 



4A;cr_ 



(^(T_ exp(Axi) + (5+ exp(— Axi))^ ' 

exp(Ax2) — 5+ exp(— Aa;2))^ ' 
and corresponding shape invariant potentials are: 

+ 



exp(Ax+) + 5^ exp(— Ax+))2 (cr_ exp(Ax_) + (5_ exp(— Ax_))^ 
^ T ^ , ^+((^-'^- exp(Ax+) - 61 exp(-Ax+)) ^ 
+ V (a_5_exp(Aa;+) + 52exp(-Ax+))2 + 

(J- exp(Aa;_) — (5- exp(— Ax_) ^ ,^^/^ 



(cr_ exp(Aa;_) + 5- exp(— Ax-))^/ \ (5_ exp(Ax2) — 5+ exp(— Aa;2))^ 

0"_ 



(34) 



((j_ exp(Aa;i) + 5+ exp(— Aa;i))2 

Explicit expressions for these potentials for particular choices of the arbitrary constant pa- 
rameters are presented in the Appendix. By the way, Eq.f lT^ for potentials f l5^ and 
has the form: 

(a + 1) — R^'^\a) = 8A^(2a + 1) ^H^^\a) + 2A^(2a^ + 2a + 1)^ . 

As for the second shape invariance ( fT71) . Eq.l fT^ allows the general solution: 

z' 

pI+p'± = A|; p± = —; 4 = ^±^±; 

(a + l)g±(a + 1) + ag±(a) = 0; g±(a) = — ; 

(X 
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with restriction for constants: 



In turn, Eq.(I7]) gives: 



+ 7-9_F = 0; (7+ + l-)f[{x,) + (7+ - 7-)/2(a^2) = 0. 



(35) 
(36) 



It is easy to check, that for 7+ 7^ ±7_ Eq. fl36|) gives /i, 2(2^1,2) ~ 2:1^25 and Eq. fl35|) can not be 
satisfied. 

Thus, the sole opportunity 7+ = ±7_ remains open. Then, either /i(xi) or f2{x2) 
vanishes, and = = A^. After all, Eq.f l35p can not be fulfilled, and no shape invariant 
potentials correspond to C± of the form (ITTl) . 
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5 Appendix: Potentials with specific parameters 
values. 

Here we give the explicit formulas of shape invariant potentials for specific values of the 
parameters. 

The general form (l32i) provides the following expressions for different choices of constants: 
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(7+ = 0, 5+ = (T_ = —S- = 1, 



^(i),(2) ^ 2AV« T 1) ^ /'exp(2Axi) + exp(2Ax2)l - h f exp(4Axi) + exp(4Ax2)(V) 
smh (Aa;_j \ / \ / 

(7+ — 0,S-^ = a- = S- = 1, 

^(,),(,) = -^^^^fef^ - h f exp(2AxO - exp(2Ax2)) - 

cosh [AX-) \ ) 

-k2 ( exp(4Axi) + exp(4Ax2) ) ; (38) 



a+ = 0+ = 



2ki + k2f 1 1 \ 2ki-k2/ 1 1 



16 \sinh (Xxi) cosh (Aa;2) / 16 \cosh (Xxi) sinh (Aa;2) 

(7+ = —5+ — a- — —5- — 1, 

^ ^ Oinh^(Ax_) ~'~ sinh^(Ax+)) 
2fci + ^2 / 1 ^ 1 \ 2A;i - A;2 / 1 ^ 1 
16 \sinh^(Axi) sinh^(Ax2) / 16 \cosh^(Aa;i) cosh^(Ax2) 

a+ — S+ — — — (5_ = 1, 



; (39) 



y(i),(2)^_2A2a(azFl)f ^ 

\cosh ( 

2ki sinh(2Ax2) + k2 



1 \ 2ki sinh(2Aa;i) + k2 



(Ax+) sinh^(Aa;_)/ 4cosh^(2Axi) 



4cosh^(2Aa;2) 



(41) 
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The analogous family of potentials for the general form (IMI) is: 

a_ = 0, 5+ = 5_ = 1, 

y(i).(2) = _4«A(2a ^ 1) |^exp(Ax+) + exp(Ax_)^ + 
+4a^|^exp(2Ax+) +exp(2Ax_)^ + 



sinh (AX2) 



cr_ = 0, 5+ = = 1, 

y(i).(2) = -4aA(2aT l)(^exp(Ax+) -exp(Ax_)^ + 

A; 



(42) 



+Aa^ I exp(2Ax+) + exp(2Ax_) ) + — — -; (43) 
\ / cosh (AX2) 



„2 



(T_ = (5_ = 5+ = 1, 

y(i),(2) k( ^ ^ \ 2aA(2a =F 1) sinh(Ax+) - 2A^a(a =F 1) + a- _^ 

\sinh^(Ax2) cosh^(Axi) / cosh^(Ax+) 
2aA(2a =F 1) sinh(Ax_) - 2A2a(a =f 1) + 



cosh^(Ax- 



(44) 



(j_ = — 5_ = 5+ = 1, 

y(i),(2) _ 2A^a(a =F 1) - 2a;A(2a =F 1) cosh(Aa;+) + 

sinh^(Ax+) 
2A^a(a =F 1) + 2aA(2a =f 1) cosh(Ax_) + 
sinh^(Ax_) 

+ 1. \ (45) 



cosh (Axi) cosh (AX2) 



cT_ = — (5_ = —5+ = 1, 

y(i),(2) _ 2A^a(a =F 1) + 2a;A(2a T 1) cosh(Ax+) + ^ 

sinh^(Ax+) 
_^2A2a(a =F 1) + 2aA(2a =F 1) cosh(Aa;_) + 
sinh^(Ax_) 

^ ' + . .1. ^ - (46) 



^sinh (AX2) sinh (Axi) 
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Some of these potentials are already known. In particular, Eq.( l37|) reproduces the two- 
dimensional generalization of Morse potential. It was investigated with essential use of its 
shape invariance in the framework of polynomial two-dimensional SUSY Quantum Mechanics 
in [12], [13], [IE]. The same conclusion concerns the potential ( l42l) after substitution a;+ = 
Vu X- = y2. 

The analogous SUSY approach was used in p!7] for the study of generalized two- 
dimensional Poschl- Teller potential, which coincides with ( 139|) after a suitable change of 
constants. 

Both potential ( HOj) and (after replacement x_ = yi, x+ = 1/2) potential (H6|) are the 
particular cases of A^— particle models of BCn type investigated in [19] in a different approach 
beyond the shape invariance context. 

To the best of our knowledge, all other cases from the above list - the potentials fl38|) . fj4T]) . 
( H3|) . fj44|) and fH5|) - are new. An additional remark is appropriate here. The transformations 
of coordinates xi,X2 to yi = X-,y2 = x+ and back link potentials ( l38l) with ( l43l) and 
potentials ( HTl) with ( l44l) . However the coupling constants in the equivalent potentials are 
different, leading to existence of double shape invariance for them, similarly to that studied 
in [15] for the Morse potential. 
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